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Q . The dynamical N = 1, SU(2) Super Yang-Mills theory is studied on the lattice 

using a new lattice fermion regulator, domain wall fermions. This formulation even at 
non-zero lattice spacing does not require fine-tuning, has improved chiral properties 
and can produce topological zero-mode phenomena. Numerical simulations of the 
full theory on lattices with the topology of a torus indicate the formation of a gluino 
' condensate which is sustained at the chiral limit. The condensate is non-zero even 

for small volume and small supersymmetry breaking mass where zero mode effects 
^ , due to gauge fields with fractional topological charge appear to play a role. 
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1 Introduction 



It is believed that super-symmetric (SUSY) field theories may play an important role in 
describing the physics beyond the Standard Model. Non-perturbative studies of these 
theories are of great interest. First-principles numerical simulations may be able to provide 
additional information and confirmation of existing analytical calculations. Typically first 
principles numerical simulations of field theories are done within the framework of the 
lattice regulator. A host of results have been produced in this way for many field theories 
and most notably QCD. Several SUSY theories can also be formulated on the lattice and be 
studied numerically. To be more specific consider the problems of putting a SUSY theory 
on the lattice (see for example [Q, ^ Q): 

1) Since space-time is discrete only a subgroup of the Poincare group survives and as a 
result SUSY is broken. This problem is not severe and is of the same nature as in QCD. The 
symmetry breaking operators that are allowed by the remaining symmetries are irrelevant. 
One can calculate at several lattice spacings a and then take the a — > limit. No fine 
tuning is needed. 

2) If the SUSY theory under consideration involves scalar fields one can have scalar mass 
terms that break SUSY since typically they are not forbidden by some symmetry. Since 
these operators are relevant fine tuning will be needed in order to cancel their contributions. 
The four-dimensional M = 1 Super Yang-Mills (SYM) theory does not involve scalars and 
therefore it does not have this problem. 

3) A naive regular izat ion of fermions results in the well known doubling problem 
For each fermion species in the four-dimensional continuum 16 are generated on the lattice 
with total chirality of zero. This results in the wrong number of degrees of freedom and 
therefore breaks SUSY. However, this problem may be possible to treat as in QCD. This 
is the case for A/" = 1 SYM. 

One possible way to remove the unwanted fermion degrees of freedom is to add an 
irrelevant operator (Wilson term ^j) that gives them heavy masses of the size of the cutoff. 
This term unavoidably breaks the chiral symmetry |Q and as a result a gluino mass term 
is no longer forbidden. Since such a term is relevant, fine tuning of the bare fermion 
mass is necessary as the continuum limit is approached in order to cancel its contribution. 
Although fine tuning is not a welcomed property this method makes it possible to recover 
the continuum target theory. 

Therefore, it is possible to simulate numerically the M = 1 SYM theory using existing 
lattice "technology" since all three difficulties can be circumvented. This observation was 
made some time ago [|l|. In particular, it was argued that, using a standard lattice gauge 
theory action with a pure gauge Wilson plaquette term and Wilson fermions in the adjoint 
representation, numerical simulations could be done. Pioneering work using these methods 
has already produced very interesting numerical results |^. Also, for proposed lattice 
tests of SYM see [^. For a supersymmetric formulation on the lattice using Kogut-Susskind 
0] fermions see [|10[] . 

There are two unwelcomed difficulties in using Wilson fermions. The first has already 
been mentioned and it is the need for fine tuning. The second is of a technical nature. It 
turns out that the Pfaffian resulting from the fermionic integration is not positive definite 
at finite lattice spacing. However, it does become positive definite as the continuum 
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limit is approached and therefore as a "cure" only the absolute value of the Pfaffian is 
used [P, 0. However, this introduces non-analyticities that may make the approach to the 
continuum limit difficult. 

Both of these difficulties can be brought under control by using an alternative fermion 
lattice regulator, domain wall fermions (DWF). The use of DWF in supersymmetric theories 
has been explored in the very nice work of @, ^. The methods in this paper are along 
the lines of these references. Domain wall fermions were introduced in were further 
developed in 



12 and in 13, 14 



They provide a new way for treating the unwelcomed 
chiral symmetry breaking that is introduced when the fermion doubler species are removed. 
Here a variant of this approach will be used [1T3| , |14|. For reviews on the subject please 
see 1 15 1 and references therein. DWF have already been used for numerical simulations 
of the two flavor dynamical Schwinger model dynamical QCD [p!7[ , 

w 
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, dynamical QCD 
as well as for simulations of 4-Fermi models 



quenched QCD 
H. The use of 



DWF in supersymmetric theories has also been explored in a different fashion in 
Furthermore, the use of overlap [1^ type fermions has been explored in 
the use of other related types of fermions has been explored in 



12, 



281, and 



In the lattice DWF formulation of a vector-like theory the fermionic fields are defined on 
a five dimensional space-time lattice using a local action. The fifth direction can be thought 
of as an extra space-time dimension or as a new internal fiavor space. The gauge fields 
are introduced in the standard way in the four dimensional space-time and are coupled 
to the extra fermion degree of freedom in a diagonal fashion. The key ingredient is that 
the boundary conditions of the Dirac operator along the fifth direction are taken to be 
free. As a result, although all fermions are heavy, two chiral, exponentially bound surface 
states appear on the boundaries (domain walls) with the plus chirality localized on one 
wall and the minus chirality on the other. The two chiralities mix only by an amount 
that is exponentially small in Lg, where Lg is the number of lattice sites along the fifth 
direction, and form a Dirac spinor that propagates in the four-dimensional space-time with 
an exponentially small mass. Therefore, the amount of chiral symmetry breaking that is 
artificially induced by the regulator can be controlled by the new parameter Lg. In the 
Lg (yo limit the chiral symmetry is exact, even at finite lattice spacing, so there is no 
need for fine-tuning. 

For the first time the approach to the chiral limit has been separated from the approach 
to the continuum limit. Furthermore, the computing requirement is linear in Lg. This is to 
be contrasted with traditional lattice fermion regulators where the chiral limit is approached 
only as the continuum limit is taken, a process that is achieved at a large computing cost. 
Specifically, because of algorithmic reasons, the computing cost to reduce the lattice spacing 
by a factor of two grows by a factor of 2^~^° in four dimensions. Therefore, the unique 
properties of DWF provide a way to bring under control the systematic chiral symmetry 
breaking effects using today's supercomputers. 

The purpose of this paper is two-fold. First, the techniques for performing a numer- 
ical simulation of the full Af = 1 SU(2) SYM theory using DWF are collected and it is 
demonstrated that they work as expected by performing numerical simulations of the full 
theory. Second, the gluino condensate is measured. It is expected that a non-zero gluino 
condensate must form ||31], |32|, |3^, Q . However, there are also arguments that the theory 
has a phase where a gluino condensate does not form P3|. In the numerical simulations 
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performed here it is found that a non-zero gluino condensate is sustained in the hmit of 
zero gluino mass. This resuh is at a finite lattice spacing and therefore SUSY is still broken 
albeit by irrelevant operators. 

It must be emphasized that in this work due to limited computer resources no attempt 
has been made to extrapolate to the continuum limit. It is possible that in this limit the 
gluino condensate may vanish. Future work using larger computer resources could calculate 
the gluino condensate at several lattice spacings and extract the continuum value. But even 
then, it will never be possible to numerically prove that the finite lattice spacing theory 
is not separated from the continuum theory by a phase transition. This problem is not 
particular to the case at hand but is in the nature of numerical investigations. They can 
provide strong evidence but not unquestionable proof. A well known case with similar 
problems relates to the question of confinement and chiral symmetry breaking in QCD. 

This paper is organized as follows. In section the DWF lattice formulation of A/" = 1 
SU(2) SYM is presented. In section 0, analytical considerations relating to the gluino mass, 
the Ward identities and the effects of topology in the patterns of chiral symmetry breaking 
are given. The numerical methods used in the simulations are discussed in section ^. The 
numerical results are presented in section ^ and the paper is concluded in section p. 

2 Lattice formulation 

In this section, the A/" = 1, SU(2) SYM lattice action and operators are presented. The 
approach is similar in spirit as to the case of Wilson fermions [jl|, ^, . The DWF formulation 
for this theory is identical to [0] and It is presented below for the convenience of the 
reader and in order to establish notation. 

The A/" = 1, SU(2) SYM theory is an SU(2) gauge theory with Majorana fermions in the 
adjoint representation. As such, the fermionic path integral results in the analytic square 
root of the corresponding Dirac determinant. This then is the Pfaffian of an antisymmetric 
matrix that has the same determinant as the Dirac operator. On the lattice, the Dirac 
operator can be defined using Wilson's approach as in |^, or the DWF approach as in 
i and i. 

The partition function is: 

Z = J [dU] J [rf^] J [d<l>]e-^ (1) 

U^{x), fj, = 1,2,3,4 is the four- dimensional gauge field in the fundamental representation, 
"^{x, s) is a (real) five- dimensional Majorana field in the adjoint representation and s) 
is a (real) five- dimensional bosonic Pauli Villars (PV) type field with the same indices as 
the Majorana field, x is the coordinate in the four-dimensional space-time box with extent 
L along each of the four directions. The boundary conditions along these directions are 
taken to be periodic for all fields. The coordinate of the fifth direction is s = 0, 1, . . . , L^ — 1, 
where Lg is the size of that direction and is taken to be an even number. The action S is 
given by: 

S = Sg{U) + Sf{^,U) + Spv{^,U) . (2) 
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Sg{U) is the pure gauge part and is defined using the standard single plaquette action 
of Wilson: 

SG = /3j:i^-lReTr[U,]) (3) 
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where /? = A/g"^ and g is the gauge coupling. 
The fermion part Spi"^, U) is given by: 

Sf = - ^{x,s)Df{x,s;x',s')-^{x',s') (4) 

x,x' ,s,s' 

where Dp is the DWF Dirac operator in the form of [0. Specifically it is: 

Dp{x, s; x', s') = 6{s - s')^x, x') + p^{s, s')6{x - x') (5) 



= W[il + ^^)V,ix)6ix + f,-x') + {l-^,)Viix')6ix' + fi-x) 
11=1 

+ (mo - 4:)5{x - x') (6) 



Pr5(1 - s') - mfPL6{Ls -1-s')- 6{0 - s') s = 

p^{s,s') = { Pr5{s + 1-s') + Pl5{s-1-s')-5{s-s') < s < - 1 (7) 

-mfPR5{0 - s') + Pl5{Ls - 2 - s') - 5{L, -1-s') s = - 1 

= '-^ (8) 

where V is the gauge field in the adjoint representation. It is related to the field in the 
fundamental representation by (see for example [Q): 

[V,ix)U = 2TT[Ulix)T'^U,{x)T'^] (9) 

and 

v,ix) = v;ix) = [v-\x)f (10) 

where T"" = ^a'^ with a'^ the Pauli matrices. In the above equations niQ is a five-dimensional 
mass representing the "height" of the domain wall and it controls the number of light flavors 
in the theory. In order to get one light species in the free theory one must set < mo < 2 



ri| . The parameter m/ explicitly mixes the two chiralities and as a result it controls the 
bare fermion mass of the four- dimensional effective theory. The dependence of the bare 
fermion mass on mo and Lg is discussed in section ^]l|. 

The fermion fleld \l/ is not independent but is related to \I' by the equivalent of the 
Majorana condition for this 5-dimensional theory [Q]: 

^=^^C'i?5 (11) 

where is a reflection operator along the flfth direction and C the charge conjugation 
operator in Eucledean space which can be set to: 

C = 7072 • (12) 
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Therefore, the fermion action can also be written as: 



J2 ^^(x,s)Mi.(a;,s;a;',s')^(x',s') 



(13) 



where 



Mp{x, s; x', s') = CR^Dpi^x, s; x', s' 



(14) 



is an antisymmetric matrix as can be easily checked 0. As a result the fermionic integral 
gives the anticipated Pfaffian: 



Because det(Ci?5) = 1 one also has that det(Mi7) = det{DF) and therefore: 



Pf(Mp 



det(L>i 



(15) 



(16) 



The Pauli-Villars action Spy is designed to cancel the contribution of the heavy fermions 



T^. Viewing the extra dimension as an internal flavor space [O] one can see that there 



are Lg — 1 heavy fermions with masses near the cutoff and one light fermion. The PV 
subtraction subtracts the Lg heavy particles. As was pointed in 0] this amounts to a 
"double" regularization of the light degree of freedom, first by the lattice and then by the 
PV field. The form of the PV subtraction used here is as in |l^ and is given by: 



Spv = "^^(^^ s)Mp[mf = l]{x, s; x', s')<l>{x', s') 



;i7) 



The integral over the PV fields results in: 



-Spv 



Pi{Mp[mf = 1]) ■ 



(1^ 



Green functions in this work are measured using four- dimensional fermion fields con- 



structed from five- dimensional fermion fields using the projection prescription [14 



X{x) = Pr^{x,0) + Pl^{x,Ls-1) 
X{x) = ^{x,L,-1)Pr + ^{x,0)Pl . 



(19) 



In the Ls oo limit of the theory these operators directly correspond to insertions in the 
overlap of appropriate creation and annihilation operators . 



Using eq. |Tl] and |l^ the Majorana condition on the four- dimensional fermion field is: 



(20) 



Because this is the correct condition for a four-dimensional field one can see that the 



definition in eq. not only produces an antisymmetric fermion matrix Mp but is also 
consistent with the projection prescription in eq. p9| as expected. 
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3 Analytical Considerations 



In this section some analytical considerations are presented. In the M = 1 SYM theory, 
a gluino mass term is the only relevant operator that can break supersymmetry and is 
also the only relevant operator that can break (at the classical level) the U{1)a symmetry. 
Therefore, the two symmetries are intimately related to the mechanisms that can introduce 
a bare gluino mass term. These mechanisms depend on the "extra" regulator parameters 
mo and Lg. This is discussed below. Next the fate of the U{1)a chiral symmetry and 
the effects of topology are presented. The chiral and supersymmetric Ward identities are 
derived in the last subsection. 



3.1 The "extra" DWF parameters 



DWF introduce two extra parameters, the size of the fifth direction Lg and the domain 
wall height or five-dimensional mass mo. These two parameters together with the explicit 



mass irif control the bare fermion mass m^s- In the free theory one finds flG 



meff = mo(2 - mo) m/ + (1 - mo) 



< mo < 2 



(21) 



In the interacting theory one would expect that mo as well as its range of values will be 
renormalized. From the above equation one can see that for the free theory the value 
of mo = 1 is optimal in the sense that finite Lg effects do not contribute to meg. In 
the interacting theory one would expect that there is no such "optimal value" since, in a 
heuristic sense, mo will fluctuate. For a more detailed analysis please see Then one 
would like Lg to be large enough so that the second term in eq. ^ will be small allowing 
for simulations at reasonably small masses and/or for dependable extrapolations to the 
mj — > 0, Ls — > cxD limit. 

The effects of finite Lg on the chiral symmetry can be best understood in the overlap 
formalism [|T2|. In that formalism a transfer matrix T along the extra direction is con- 
structed. Because the gauge fields are not changing along that direction the product of 
transfer matrices simply results in T^". For Lg = oo this is a projection operator that 
projects the reference vacuum state to a ground state. The fermion determinant is then 
the overlap of the reference vacuum state with that ground state. In |r2| it was shown that, 
as a lattice gauge field configuration changes, from say the zero topological sector to sector 
one, an eigenvalue (or a degenerate set of eigenvalues) of the corresponding Hamiltonian 
H changes sign. As a result, the filling level of the ground state becomes different from 
that of the reference vacuum state. Then the overlap is zero indicating the presence of an 
exact zero mode. This remarkable property is maintained to a good degree even at finite 
Lg as was found in [^. Unfortunately, this property is also the reason for most of the 
difficulties with DWF. As the eigenvalue of the Hamiltonian H changes sign it crosses zero. 
In such a configuration the transfer matrix has an eigenvalue equal to one and therefore 
even at = oo there is no decay along the extra direction, the two chiralities do not 
decouple, and chiral symmetry can not be restored. Fortunately, configurations for which 



H has an exact zero eigenvalue (for a given mo) are of measure zero |jT2|, and therefore 
are of no consequence. However, configurations in their neighborhood are not of measure 
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zero and such configurations will exhibit very slow decay rates. Therefore, in order to 
restore chiral symmetry, very large values of Lg may be needed. Since one would expect 
that the neighborhoods of such configurations are suppressed closer to the continuum limit 
this problem should become less severe as that limit is taken. This has been observed in 
numerical simulations of the Schwinger model ||T6|, of full QCD [O], and of quenched QCD 



m 19, 22, 23 



In the region where it makes sense to parameterize these effects by a residual mass in 
an effective action it has been found that: 

mcff = Com/ + m^es, m^s = Ci exp(-C2Lj (22) 

where for dynamical QCD at the currently accessible lattice spacings the decay is found to 



be C2 ~ 0.02 IjT^. For quenched QCD the situation is better because current computing 



resources can simulate lattices with smaller lattice spacing. There, a value of C2 ~ 0.1 is 



found |T8|, |T9|, Also in these studies the value of C2 was a weakly changing function of 
mo indicating that for practical purposes there is no optimal value of mo. 

In the case of the M = 1 SYM SU(2) theory the Hamiltonian corresponding to the 
five dimensional transfer matrix has eigenvalues that are doubly degenerate because the 
fermion fields are in the adjoint representation Therefore when there is a "topology" 



change two eigenvalues will have to cross through zero (as compared to one for fundamental 
fermions). This may make this theory harder to study than QCD in the sense that larger 
Ls values may be required. On the other hand, since no massless Goldstone particles are 
expected, the sensitivity of the spectrum on Lg may be considerably milder. In any case, 
in this paper the only fermionic observable that will be discussed is the gluino condensate. 
This quantity is known to approach its Lg —>■ oo limit with faster decay rates than the ones 
in rrij-es (for a discussion and results for full QCD see [|17|; there the decay rate for the chiral 
condensate was about five times faster than that for mres)- 

As was discussed above, the range of mo is renormalized by the interactions. It has 
been found that as the lattice spacing increases and one moves away from the continuum 
limit this range shrinks in size and for currently accessible spacings in QCD that range is 
about [1.4,2.0]. As one moves even farther away from the continuum limit this range can 
shrink to zero and then it will not be possible to have light DWF modes [0, However, 



it must be emphasized that for as long as the range of allowed values of mo is not of zero 
size the overlap formalism, although it does not specify how it is approached, guarantees 
the existence of the Lg —>■ oo limit. In this work, mo = 1.9 and, as it will be shown in 
section 5, the behavior of the gluino condensate vs. Lg is consistent with an exponential 
ansatz. 

3.2 Chiral symmetry and topology 

Fermions in the adjoint representation of the SU(N) gauge group have a Dirac operator 
with index: 

2Ni^ (23) 

where u is the winding of the background field configuration. Classical instantons have 
integer winding and they cause condensation of operators with 2A^ Majorana fermions. 
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This results in the breaking of the U{1)a chiral symmetry down to the Z2N symmetry 
by the corresponding anomaly. The remaining Z2N symmetry may break spontaneously 
down to Z2 [0. Mechanisms for this further breaking have been explored for example in 
where instantons and fractionally charged objects such as torons or caloron 



monopole constituents [^] were investigated as the source of this symmetry breaking. 

Since in a toroidal geometry fractional winding numbers are possible |37|, the partition 
function of the full theory can be expressed as: 

^W = Ee"%> z/ = 0,±l/iV,±2/iV,... (24) 

where Q is the vacuum angle and is the partition function on the sector with winding v. 
For the theory with a soft breaking by a mass m/ the interplay of the volume and mass in 
the formation of the gluino condensate has been analyzed in The reader is referred to 
that reference for a very nice presentation on the subject. Assuming a mass gap is present 
in the theory the authors of |^0[ show that non-zero contributions to the gluino condensate 
(xx) come almost exclusively from the v = 1/N sector H ntf 'xV x (xx) <^ 1- On the other 
hand, ii rrif xV x (xx) ^ 1 ^ill sectors contribute to a non-zero condensate. 

The above considerations result in an unusual picture. If the infinite volume limit is 
taken (followed by the massless limit) it is possible that a gluino condensate will form due 
to spontaneous breaking of the discrete symmetry Z2N down to Z2. On the other hand, at a 
finite volume and zero mass a gluino condensate can form due to the presence of fractional 
winding configurations. Since the volume is finite, this can not be the result of spontaneous 
symmetry breaking. Instead, it is similar to symmetry breaking due to topological effects 
as, for example, in one flavor QCD. As pointed above the size of mj x V" x (xx) controls 
which "scenario" takes place. 

On the lattice there is no clear definition of topology. The path integral over the SU(N) 
group space generates configurations of all possible windings. In order for the lattice theory 
to be able to reproduce phenomena that relate to topology it is essential that the lattice 
Dirac operator obeys the index theorem in a statistical sense. This is highly non-trivial 
since it is obviously related to the doubling problem. Traditional fermions (Wilson or 
staggered) do not exhibit exact zero modes at finite lattice spacing. On the other hand. 



as mentioned in section 3.1, DWF at L,, = 00 have exact zero modes and at finite 



have robust zero modes to a good approximation PD[. An approximate form of the index 



theorem has been found to be obeyed for fundamental fermions in the overlap formulation 
in quenched SU(2) |T| and in quenched SU(3) |g. 

The index of adjoint fermions in the overlap formulation in quenched SU(2) has been 
studied in |^^. In that work it was pointed out that the overlap Dirac operator for adjoint 
fermions in the SU(2) gauge group is necessarily even- valued. Then the question posed by 
the authors of pSf was whether or not all even values are realized or only values that are 
multiples of four are present. The latter case corresponds to configurations with instantons. 
The former case corresponds to fractional winding numbers. Configurations with fractional 
winding were found and their presence persisted as the lattice spacing was decreased 

In this paper DWF are used at finite Lg and therefore some of the clarity present in the 
Lg = 00 case will be lost. However, the full theory (including the fermion determinant) is 
studied here. Furthermore, it is interesting to see if at a small volume and zero mass the 
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gluino condensate still forms and if it does to what extent its value is due to zero mode 
effects. The numerical results are presented in section 



3.3 Ward identities 

As discussed in the introduction and in section the DWF formulation of the A/" = 1 
SU(2) SYM theory at the Lg ^ oo limit is expected to preserve the U{1)a chiral symmetry 
(at the classical level) and break supersymmetry only by irrelevant operators. Since the 
DWF formulation contains many more fields than the continuum theory, one may naturally 
wonder what are the SUSY transformations in terms of these fields. In particular, while 
the continuum theory has a single Majorana fermion the DWF lattice theory contains Lg 
Majorana fermions and Lg corresponding PV fields. Since all these fields, except for one 
Majorana fermion, have masses near the cutoff, one can expect that the SUSY transfor- 
mations should only transform the gauge field and the light Majorana fermion represented 
by the boundary field x of sq. 0. Similarly, the chiral symmetry transformations should 
only involve the field x- However, one should expect that this choice of SUSY and chiral 



transformations is not unique. For example, see for a different choice of QCD chiral 
transformations that involve all fermion fields in one half of the fifth direction transforming 
vectorially and all fermions in the other half also transforming vectorially but with opposite 
charge. That choice could also be appropriate here for the chiral transformations, but it 
may make the SUSY ones more complicated. 

As a first step in deriving the Ward identities, the fermionic part of the action in eq. § 
is rewritten in terms of the boundary field x'- 

Sp = + Sp^ (25) 

where Spo does not depend on the field x 

Sp^ = - i xix)pNix, x')xix') - mfxix)6{x, x')xix') 

—x{x)B{x,x')(f){x') — (f){x)B{x,x')x{x') (26) 

where 

0(x) = PR^{x,Lg-l) + PL^{x,0) 

0(x) = ^ix,0)Pn+^ix,Lg-l)PL (27) 

are the "wrong" projected fields in the sense that they are defined on the opposite wall 
from where the corresponding light mode is localized. If indeed there is localization one 
would expect that in the Lg —>■ oo limit these fields will have no overlap with the light 
mode. The operator is the naive part of the four-dimensional Wilson operator in eq. ^ 



and B is the symmetry breaking part {B is the equivalent of B in [|12| , |T4|): 



1 ' 



Pn{x,x') = - X! 7m V^{x)S{x + fi- x') - V^{x')S{x' + fi - x) (28) 



10 



1 ^ 

B{x, x') = (5 - mo)S{x -xi)--Y, \y^{x)6{x + /i - x') + V^^{x')6{x' + /i - x)] . (29) 



These operators have the following properties: 



0, 
0, 



= B^ = B 



(30) 



First the Ward identity corresponding to the U{1)a symmetry is derived. The symmetry 
transformations are: 



ia{x)'j5x{x) 
ia{x)xix)-f5 



(31) 



where a{x) is an infinitesimal real number and 6a symbolizes the change under the chiral 
transformation. Then the Ward identity is: 

{AMx)Oiy)) = 2mf{Ux)0{y)) + 2{JB{x)0{y)) + i{5AO{y)) (32) 

where the backward difference is defined as A^/(a;) = /(x) — /(x — /x). The currents are: 

J^,{.x) = xix)-f5l^V^ix)xix + /i) (33) 

M^) = xix)-f5Xix) (34) 

Jb{x) = J2x{xh5B{x,y)(i){y) . (35) 
y 

If in the above Ward identity 0{y) = J^iy) one gets 

/^,{J,{x)Uy)) = 2mf{Mx)My)) + 2{JB{x)My)) " 2(x(y)x(y)) (36) 

In this identity the term with Jb will be responsible for producing the ABJ anomaly in 
the Ls — > cxD limit. On the other hand, if Lg is kept finite this term is similar to the 
one for Wilson fermions which, besides producing the ABJ anomaly, also produces a mass 
redefinition. For an analysis of QCD with DWF at finite Lg see [p!9[ . 

As mentioned earlier these chiral transformations are different than the ones in [Q. If 
the transformations relevant for a non-singlet current in QCD were done on the fields x^ 
X, one obtains a Ward identity exactly as in |lj] but with the currents A°(x) and J^g{x) 
replaced with: 



A'.ix) 



x{x)-f5l,.^^U^{x)xix + + xix + /x)757/.A"[/^(x)x(a:) 



(37) 



^ [xix)^5X^B{x,y)<P{y) + <P{yh,X^B{y, 



x)X[x) 



{3t 



The derivation of the SUSY Ward identity is similar to the one for Wilson fermions. One 
can use the existing calculations for Wilson fermions |1], 0, ^ to elucidate the differences 
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between the two formalisms. Here the derivation in [H3] will be followed. The symmetry 



transformations are as in |^3[ and commute with parity. 

The change of the pure gauge action with respect to the transformation of the gauge 
field is of course the same. In terms of the symmetry breaking part of the Ward identity it 



contributes a term denoted below by X2{x) +X3{x) where X2, X3 are as in This term 



breaks SUSY because of the explicit breaking of the Lorentz symmetry. Using improved 
pure gauge lattice actions can alleviate the effects of this breaking. Such an improvement 
is not considered here. 

The change of the fermion and Pauli-Villars parts of the DWF action with respect to 
the transformation of the gauge fields produces terms for all Lg slices. In particular the 
variation of the fermion matrix Dp of eq. ^ with respect to the gauge field is: 

6uDp{x,x'; s, sf){mf) = 6{s — s')6up{x,x') (39) 

One sees that 6Dp is independent of mj and is diagonal in the fifth direction. Furthermore 
S]p{x, x') is the same as the variation of the Wilson operator. Therefore, this variation 
contributes to the symmetry breaking part of the Ward identity the terms: 

Xf(x) = VXf(x,s) (40) 



and 



Xr(a;) = Vxr(x,.) (41) 



where X^ (x, s) is as X4 in except that the four- dimensional Wilson fermion fields that 



have their spin indices contracted are replaced by \I^(x, s), \l/(x, s) while the other Wilson 
fermion field is replaced by x(^)- Similarly X^^(x, s) is as X4 in |4^ except that the four- 
dimensional Wilson fermion fields that have their spin indices contracted are replaced by 
the Pauli-Villars fields $"^(x, s)CR^, $(a;, s), the other Wilson fermion field is replaced by 
x{x) and the sign of the second term in X4 is minus instead of plus due to the commutativity 
of the Pauli-Villars fields. 

The change of the action with respect to the fermion field transformations can be par- 
tially deduced from the corresponding Wilson fermion calculation. Since this transforma- 



tion only involves the action Sp^ in eq. ^ one can observe that the first two terms of that 
action are identical with the action of naive fermions (Wilson fermions with r = 0). These 
will contribute identical terms as the r = part of the Wilson action. They contribute 
to the divergence of the SUSY current and to the mass term of the Ward identity given 
below. Finally, the transformation of the last term of the action Sp^ in eq. ^ is easy to 
calculate and is denoted by Xi{x): 

X,{x) = -P;^{x)ap^B''\x,x')<p\x') . (42) 



This term is closely related to Xi of |^ 
The Ward identity is: 



{A,S,{x)0{y)) = m;{Dg{x)0{y)) + {Xs{x)0{y)) - {6sO{y)) (43) 
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where the supersymmetric current and the quantity Ds are as in [43|. The symmetry 



breaking term Xs{x) is also similar to the one in 



Xsix) = X,{x) + X2ix) + Xsix) + Xi (x) - ^ (x) (44) 

As mentioned above the symmetry breaking term X2{x) + X^{x) is due to the breaking 
of Lorentz symmetry by the lattice. The X^{x) and X^^[x) terms break the symmetry as 
in Wilson fermions. These terms do not cancel each other exactly Q However, one would 
expect large cancellations of heavy modes. The terms in Xl[x,s) that are proportional 
to the Wilson parameter involve fields that couple to the light modes by an amount that 
is exponentially small in Lg. One would expect these terms to be nearly canceled by 
the corresponding Pauli-Villars terms resulting in exponentially small contributions. The 
remaining terms that involve fields away from the relevant domain walls should also yield 
similar cancellations. As a result the only terms that should make significant contributions 
should be the ones that involve fields of the "correct" chirality near the domain walls. These 
few terms would couple to the light modes and be further regularized by the corresponding 
Pauli-Villars terms. Clearly this analysis of cancellations is heuristic. A detailed calculation 
using for example perturbation theory or transfer matrix methods would be interesting but 
it is beyond the scope of this paper. 

Finally, the symmetry breaking term Xi{x) involves the field (t){x) that is expected 
to have no overlap with the light mode in the Lg ^ oo limit. If Lg is finite then DWF 
are similar to Wilson fermions and an analysis as in should indicate that this term is 
responsible for the same mass redefinition as the one in the chiral Ward identity. 



4 Numerical methods 

As can be seen from section |^ the A/" = 1 SU(2) SYM theory can be simulated as a theory 
with 0.5 flavors of Dirac fermions in the adjoint representation. An efficient and popular 
algorithm that can be used to simulate any number of flavors is the hybrid molecular 



dynamics R (HMDR) algorithm of |^^. Because of the Grassmann nature of fermions 
these algorithms need to invert the matrix Dp of eq. |. That matrix is not Hermitian. 
This is a problem since some of the more efficient inversion algorithms require the matrix 
to be Hermitian. However, because: 

^^RDpRi^ = D^F (45) 

one has that: 

detlDpf = detlDpDl,] . (46) 

Then one can invert the Hermitian matrix DpD^p and then use the HMDR algorithm to 
take the appropriate power so that the desired number of flavors is simulated. This method 
is adopted here and the 0.25 power is taken in order to go from a theory with two Dirac 
fermions to a theory with one Majorana fermion. In other words, the fermion determinant 
that is used in the simulation is: 

(det[DirD^])°-25 = I det[Dir]|°-^ = dei[Dpf '' (47) 



^We thank Y. Shamir for pointing this out to us 
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where in the last equahty use was made of the fact that for non-negative rrif detlDp] is 
also non-negative 0. This approach was also taken in for Wilson fermions. For an 
approach that uses Wilson fermions and the multibosonic algorithm instead of the 
HMDR algorithm see . However, as mentioned earlier in the case of Wilson fermions the 
last equality in eq. is not true for all gauge field configurations. 

The HMDR algorithm uses molecular dynamics methods in order to produce the correct 
statistical ensembles. Because the molecular dynamics step size 6t is finite discretization 
errors are introduced. There are two ways one can deal with this problem. One is to 
simulate at various values of 6t and then extrapolate to 6t = 0. Another method is to use 
St small enough so that the errors are negligible when compared with the statistical errors. 

In order to ensure this, one can simulate the two Dirac flavor theory at the same 
parameters and same 6t. For the two flavor theory, one has a local action and therefore, at 
the end of the evolution, one can employ a Metropolis accept-reject step. Then the finite 
6t errors are "converted" to a non-ideal acceptance rate and in effect they are reflected in 
the final statistical errors. This is the exact hybrid Monte Carlo $ (HMC$) algorithm of 
PB| , P[ . Therefore the acceptance rate is an indication of the size of the finite 6t errors 
in the HMD integration. By simulating the two Dirac flavor theory with (HMC$) one can 
set 6t so that the acceptance rate is high, say ~ 90%. Since the coefficient of the finite 6t 
errors is proportional to the number of flavors one would expect that for 0.25 flavors the 
errors would be small and at the few percent level. 

The only fermion observable measured in this work is the gluino condensate. By in- 
serting appropriate source terms as in the gluino condensate was measured as the trace 
of Dp^ with spin and fifth-direction indices restricted as dictated by eq. |TP| The trace was 
calculated using a standard stochastic method. All inversions in this work were done using 
the conjugate gradient (CG) algorithm. An even-odd preconditioned form of the matrix 
was inverted. For more details on the numerical algorithms and methods employed 
to DWF simulations see |jl6|, |17|. 



5 Numerical results 
5.1 Simulation parameters 

In all simulations the domain wall height was chosen to be mo = 1.9. As mentioned in the 
previous section, the finite St errors were kept to the few percent level by using a small 
6t. For all simulations the step size was set to 6t = 0.01 and the trajectory length to 
r = 0.5. In order to confirm that this choice introduces finite step size errors that are small 
compared to the statistical errors an HMC$ simulation for two Dirac flavors was run for 
Ls = 12 and ruf = 0.04. It produced an acceptance rate of ~ 90% suggesting that the 
finite 6t errors of the 0.5 flavor theory are small. Furthermore, an HMDR simulation was 
also run for two Dirac flavors using the exact same parameters. The value of the gluino 
condensate obtained from these two simulations was the same within statistical errors. 

The CG stopping condition for all simulations was set to 10~^ for the evolution and 
to 10~^ for the calculation of XX- The number of CG iterations varied between ^ 100 for 
TUf = 0.08, Ls = 12 and 250 for ruf = 0.0, = 24. 
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The 8^ volume simulations were done with (3 = 2.3. The value of (3 was chosen so that 
one is not close to the point where the box size becomes too small and a thermal transition 
takes place, but also not too deep in the strong coupling regime where the finite effects 
become severe. The transition point of the Nt = 8 quenched theory is at /5 = 2.5115(40) 
P7|| . In figure 1 the magnitude of the fundamental Wilson line measured in quenched 

simulations in an 8^ volume is plotted vs. (3. In the quenched theory this is an order 
parameter. As can be seen from that figure, a rapid crossover takes place around jS = 2.5. 
In the same figure the value of {\W\) from a simulation of the dynamical theory at /3 = 2.3 is 
also shown (cross). The quenched and dynamical values are very similar indicating that at 
f3 = 2.3 the dynamical theory is in a phase that "confines" fundamental sources. Therefore, 
the box size is large enough to avoid finite temperature effects that would of course spoil 
SUSY. Using the quenched theory as a guide the 4*^ simulations were done at /3 = 2.1 since 
the quenched transition at A^^ = 4 is known to take place at /? = 2.2986(6) |^^. At /? = 2.1 
the lattice spacing is larger than at (3 = 2.3. However, the lattice sizes are small and do 
not allow a reliable measurement of the lattice spacing. According to (3 = 2.1 — 2.3 is 
in the beginning of the weak coupling regime. Then if one uses the weak coupling form in 
48| one finds that the lattice spacing at /5 = 2.1 is about a factor of two larger than the 



one at /5 = 2.3. 

In order to estimate the necessary number of thermalization sweeps two simulations 
were run on an 8^ lattice at f3 = 2.3, Ls = 12 and mf = 0.04. The first simulation used an 
initial configuration with all gauge links set to the identity (ordered) and the other used an 
initially random configuration (dis-ordered). The evolutions in "computer time" are shown 
in figure 2. As can be seen, the two ensembles converged after about 100 sweeps. This 
number of thermalization sweeps was then used in all other simulations which were started 
from an ordered initial configuration. The number of measurements after thermalization 
for all simulations is about 200 with measurements done in every trajectory. The gluino 
condensate was measured with a single "hit" stochastic estimator. 

5.2 The gluino condensate at the chiral Hmit 

In order to be able to extrapolate to the chiral limit, corresponding to — >■ oo and mj = 0, 
the mass ruf and the size of the fifth direction Lg was varied. The results of all simulations 
are given in tables |l| and § in the appendix. Three different methods were used to analyze 
the data and calculate the gluino condensate in the chiral limit. 

I. For fixed Ls, the data for m/ = 0.08, 0.06, 0.04, 0.02 were fit to a function: 

bo + birrif . (48) 

This functional form is valid provided m/ is small enough. Otherwise, higher order terms 
must also be included. The data and fits are shown in figure 3 and the results of the fits 
are given in table ^. Then the extrapolated values bg were fit to a form: 

Co + ciexp{-C2Ls) . (49) 

This functional form is approximate but it is expected to be valid close enough to the 



continuum and has been found to be consistent in simulations of the Schwinger model [16 
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and of QCD even at relatively large lattice spacings (see for example |17|). The data and 
fit is shown in figure 4 and the results of the fit are given in table g 

II. For fixed ruf the data for Lg = 12, 16,20,24 were fit to the form of eq. The 
data and fits are shown in figure 5 and the results of the fit are given in table H. Then the 



extrapolated values cq were fit to the form of eq. |^. The data and fit is shown in figure 6 
and the results of the fit are in table ^. 

III. Additional simulations were done for rrif = and L<, = 12, 16, 20, 24. The data were 



fit to the form of eq. ^ The data and fits are shown in figure 7 and the results of the fit 
are in table |^. 

The — i> and oo extrapolated values of the gluino condensate for each one 

of the above three methods are summarized in table ^ As can be seen, all values are 
consistent within the statistical errors. This suggests that the systematic errors inherent 
to the limited statistics and to fits onto functions that represent the data only for a limited 
range are small. Furthermore, it suggests that the fitting functions used are consistent 



(please see subsection ^.4| for more discussion on the validity of these fitting functions) 



5.3 The telltale signals of topology in numerical simulations 



In order to investigate the issues discussed in section the gluino condensate was also 
calculated in a smaller 4^ lattice volume at /? = 2.1. It was measured only for rrif = and 
method III above was used to extrapolate to the — > oo limit. The data and fit are shown 
in figure 8 and the fit results are given in table ^. The 8^ data from figure 7 are presented 
again in this figure to aid comparison. The value has decreased indicating that scaling is 
violated. However, without more simulations at other lattice spacings and volumes one can 
not conclude much from this result. The f3 = 2.1 coupling is in the strong coupling region 
and furthermore the 4^ lattice volume is rather small. 

However, it is interesting to notice that the parameter V x {xx) l^^oo ~ 8.4 (a factor of 
12 coming from the normalization of (xx) has been included). Since m/ = the effective 
mass meff gets its value from finite Ls effects. As Ls is increased rries becomes small. From 
analysis of rrics in strong coupling QCD [1^ one would roughly guess that rries < 0.1. Then 



[meff X X (xx) Ls^oo] < 1- In that case, the analysis of |^ can be followed and one would 
expect the value of the condensate in the 4"^ lattice to be mostly supported by configurations 
with total winding of ±1/2. Indeed, this can be seen from figure 9. In that figure the 
evolutions in "computer time" are shown. The "spikes" in the evolution are apparent and 
they become more pronounced and less frequent as Lg is increased (and in effect rries is 
decreased). This is exactly how the effect of zero modes for small [rrtefj x V x (xx)ls-^oo] 
would present itself in a numerical simulation of the dynamical theory. As the fermion mass 
is made smaller, (xx) is expected to receive most of its value from sectors with winding 
±1/2. However, in these sectors the fermion determinant is very small because of the zero 
mode. Since the probability for the algorithm to generate a gauge field configuration is 
proportional to the fermion determinant one would expect that these sectors will be visited 
less and less frequently as the effective mass is decreased. When these sectors are visited the 
value of XX will be very large (spikes) in order to compensate for the infrequent sampling. 
In this way the presence of the zero mode in the observable "balances" the presence of the 
zero mode in the determinant. As the mass is made smaller one would have to increase 
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the size of the statistical sample in order to include enough of these increasingly "rare" but 



very large fluctuations. For similar results in the Schwinger model and QCD see |T^, . 

A histogram of the values of XX is presented in figure 10 (solid line). For small Lg the 
effective mass is larger and xx is distributed with a symmetric looking distribution around 
the mean value. However, for Lg = 40 the effective mass is smaller and the distribution 
has a more pronounced "tail" towards larger values. In order to investigate this further 
numerical simulations at exactly the same parameters, but without the fermion determi- 
nant (quenched theory) were done. The histograms from these simulations are shown in 
the same figure for comparison (dotted lines). One can observe that the absence of the 
fermion determinant had the effect of shifting the distributions to higher values. This is 
expected since configurations with small eigenvalues that produce larger values of xX ci-re 
not suppressed anymore and are produced more frequently. Also, one can observe that the 
number of configurations with xX larger than ^ 0.007 that appeared as spikes in figure 9 
have now increased in number. These observations lend support to the presence of small 
near-zero eigenvalues. Furthermore, configurations with fractional topological charge have 
already been found in quenched SU(2) simulations at similar couplings [^. It would be 
very interesting to calculate the index for the configurations of figure 9 using the methods of 
and see to what extent there is a correlation between fractional topological charge and 
the observed spikes. This correlation should be exact at ^ oo but it will be obscured 
at finite Lg by the presence of non-zero nies. This investigation is beyond the scope of this 
work. 

Furthermore, it should also be noted that on the 8'^ lattice there are no visible spikes 
up to Lg = 24. This can be seen in figure 11. Presumably this is because the product 
[rries xV X {xx) Ls^oo] is probably much larger than in the 4^ lattice. Again, this statement 
is not exact since the value of meg was not measured. 

These results are consistent with the discussion in section p.2| . However, since even 
with rrif = 0, sua Lg extrapolation is essentially an extrapolation from non zero masses 
these results are not necessarily the results of a simulation at exactly Lg = oo. It is 
still possible that if such a simulation were done one could have found that the gluino 
condensate is zero. This could happen since in a finite volume and zero mass the effects of 
spontaneous symmetry breaking are absent and the zero mode effects alluded to above may 
not be sufficient to sustain a non zero vacuum expectation value. This type of simulation 
is possible and can be done using the overlap formalism [12] or exact Neuberger fermions 
0]. However, if one is to maintain exact chiral symmetry these methods will demand large 
computing resources. 



5.4 The fine print 

Perhaps the largest uncertainty in the analysis presented in the previous subsections has to 
do with the assumption of exponential decay as in eq. ^ For small enough lattice spacings 
and large enough Lg this behavior is expected to be true. All data presented in this work 
were well represented by this ansatz. However, as with any numerical investigation, one can 
never completely disprove all other possibilities. While such an exercise over all possible 
functions would clearly be fruitless there are some alternative forms that may be reasonable 
to consider since they are based on analytical considerations. 
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Far from the continuum hmit, the approach to the chiral hmit may become power law 
4^ or even completely disappear p|, |39|. In order to explore the possibility of power law 



behavior the ruf = data for the 8^ and 4^ volumes were fit to the form: 

do + diLf . (50) 

The results of the fit are given in table ^ (the fits are not presented in any of the figures). 
As can be seen from that table the x^/dof of these fits is significantly larger than the one 
of the corresponding exponential fits to the same data. 

Another possibility is decay to zero with two different exponential decay rates. Such a 
behavior was found to be consistent with investigations of the two flavor Schwinger model 
p!6| for a quantity that is expected to vanish in the chiral limit. There it was argued 
that the fast decay rare is due to fluctuations within a given topological sector while the 
slow decay rate is due to the presence of topology changing configurations. Therefore, for 

= one could try to fit the largest three points to the form: 

eoexp(-eiLj . (51) 

The results of the fit for the 8^, j3 = 2.3, mj = 0, and Lg = 16, 20, 24 points as well as for 
the 4"^, P = 2.1, nif = 0, and = 24, 32, 40 points are shown in table 0. The fit to the 8^ 
data is acceptable. However, the fit to the 4^ data has a rather large x^/dof. Because this 
fit is for larger Lg than the 8'^ fit one would expect that if there were a second exponential 
decaying to zero its effect would be more pronounced in the 4^ fit. Therefore, the large 
X^/dof of the 4"^ fit suggests that the presence of a second exponential decaying to zero is 
not likely. This could be made more precise if simulations with larger Ls values for the 8^ 
and 4^ lattices were done. However, such simulations are beyond the computing resources 
of this project. Also, the analysis in suggests functional forms with more parameters. 



It would be interesting to fit to these forms but that would require more data points and 
better statistics both of which are also beyond the computing resources of this project. 

Finally, the SUSY breaking by the irrelevant terms may have non-negligible effects at 
the lattice spacings studied here. Although it was found that the chiral condensate is non- 
zero at the chiral limit in two lattice spacings, this is not enough to estimate its value in 
the continuum limit. 



6 Conclusions 

The formulation of A/" = 1, SU(2) Super Yang-Mills theory on the lattice with domain wall 
fermions (DWF) has several advantages over more traditional lattice fermion regulators. 
Even at non-zero lattice spacing the chiral limit can be taken by letting Lg oo, where 
Lg is the number of sites along the fifth auxiliary direction. Since in that limit there is no 
gluino mass term, supersymmetry is broken only by irrelevant operators and there is no 
need for fine tuning. Also, in that limit the theory has exact zero modes on non-trivial 
topological backgrounds. 

However, even at finite Lg, where numerical simulations are done, these properties are 
maintained to a good degree allowing extrapolations to the Lg oo limit. Furthermore, 
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the Pfaffian resulting from the integration of Majorana fermions is positive definite at finite 
Lsi non-zero lattice spacing and for any background gauge field configuration. As a result, 
one can unambiguously interpret it as a probability measure to be used by the numerical 
simulation for importance sampling. This property also allows the use of standard numerical 
algorithms where any number of flavors Nf can be simulated. By contrast, Wilson fermions 
have this positivity property only at the continuum limit. 

In this work, the full A/" = 1, SU(2) Super Yang-Mills theory was numerically simulated 
on the lattice using DWF. The gluino condensate (xx) was measured. These simulations 
did not present any unexpected technical difficulties. 

A flnitc vahie of Lg breaks chiral symmetry and induces a small gluino mass. In addition, 
an explicit gluino mass rrif was used to provide extra control. Several and Lg values were 
used (all corresponding to positive gluino mass) and the value of (xx) was extrapolated to 
the chiral limit using three different methods. All methods gave consistent results indicating 
small systematic effects and suggesting that the functions used for the fits are consistent. 
These simulations were done on a lattice with 8^ lattice sites. 

Additional simulations on a lattice with 4^ lattice sites but approximately double the 
lattice spacing were done. Again, extrapolations to the chiral limit gave a non-zero (xx)- 
In this lattice [mass x volume x (xx)] < 1- Then analytical considerations suggest that 
the value of (xx) must come mostly from topological sectors with fractional topological 
charge of ±1/2. Indeed, as the mass was made smaller unusually large values (spikes) 
were observed in the statistical sample of xx indicating the singular contribution of these 
sectors. 

The spectrum of the theory is of great interest but it was not possible to measure on 
the small lattices considered here. Also, the gluino condensate was measured only on two 
different lattice spacings and therefore it was not possible to extrapolate to the continuum 
limit where comparisons with analytical results would be possible. Future work could 
explore these very interesting topics. 
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7 Appendix 

In this appendix the various tables are presented. 
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Table 1: The values of (xx) for the 8^ simulations a,t P — 2.3, mo = 1.9. 
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Table 2: The values of (xx) for the 4"^ simulations at /3 = 2.1, mo = 1.9. 
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Table 3: The results of the linear fits presented in the various figures to the function 
bo + biuif (dof is short for degree of freedom). 
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Table 4: The results of the exponential fits presented in the various figures to the function 
Co + Cl exp(— C2-Ls) (dof is short for degree of freedom). 
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{xx)imf ^ 0,L, ^ oo) 
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II 
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Table 5: The m/ — > 0, — > oo extrapolated values of (xx) using the three different 
extrapolation methods described in the text. 
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Table 6: The results of the power law fits in the data of figures 7 and 8 to the function 
do + diLg"^ (dof is short for degree of freedom). 
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Lattice 
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Table 7: The results of the exponential fits without a constant using the function 
eoexp(— eiLs) (dof is short for degree of freedom). 
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Figure 1. The magnitude of the Wilson hne < \W\ > in an 8^ lattice. The diamonds 
are from the quenched theory. The cross is from the dynamical theory ai P — 2.3 with 
Lg = 24, ruf = 0.0 and mo = 1.9. 
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Figure 2. The gluino condensate values generated by the computer starting from an 
ordered initial configuration (solid line) and from a dis-ordered initial configuration (dotted 
line). The x axis is the configuration number and corresponds to "computer time". This 
is from a simulation of the full theory on an 8^ lattice at /5 = 2.3 with = 12, ruf = 0.04 
and mo = 1.9. 
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Figure 3. The gluino condensate vs. rrif for various values of Lg from tlie dynamical 
theory on an 8^ lattice aX (3 = 2.3 and mo = 1.9. The fits are to the function in eq. 
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Figure 4. The extrapolated values of the fits in fig. 3 vs. L^. The fit is to the function in 
eq. |9. 
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Figure 5. The gluino condensate vs. Lg for various values of mj from tlie dynamical 
theory on an 8^ lattice aX (3 = 2.3 and mo = 1.9. The fits are to the function in eq. 
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Figure 7. The gluino condensate vs. Ls from the dynamical theory on an 8^ lattice at 
j3 = 2.3, frif = 0.0 and rriQ = 1.9. The fit is to the function in eq. 
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Figure 8. The gluino condensate vs. from the dynamical theory on a 4^ lattice at 
j3 = 2.1, ruf = 0.0 and ttiq = 1.9 (crosses). The 8^ data of fig. 7 are also plotted for 



comparison (diamonds). The fits are to the function in eq. ^ 
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Figure 9. The gluino condensate values generated by the computer for the full theory on 
a 4^ lattice at f3 = 2.1 with rrif = 0.0, mo = 1.9, and for various values of Lg. The x axis 
is the configuration number and corresponds to "computer time" . 
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Figure 10. The histogram of the values of the gluino condensate from fig. 9 (sohd 
hnes). The dotted hnes are from an identical set of simulations, except that the fermion 
determinant has been set to one (quenched theory). The time evolutions for the quenched 
simulations are not shown in this paper. The area under the curves is normalized to one. 
All sets have 300 data points. 
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Figure 11. The gluino condensate values generated by the computer for the full theory 
on an 8^ lattice at P — 2.3 with mj = 0.0, mo = 1.9 and Lg — 24. The x axis is the 
configuration number and corresponds to "computer time" . 
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